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ABSTRACT 

We present an evolutionary picture of a forming star. We assu me a singul ar, isothermal sphere as the 
initial state of the core that undergoes collapse as described bv iShiJ dUza). W e include the evolutio n 
of a first hydrostatic core at early times and allow a disk to grow as predicted bv I Adams & 

We use a 1-dimensional radiative transfer code to calculate the spectral energy distribution for the 
evolving protostar from the beginning of collapse to the point when all envelope material has accreted 
onto the star-bdisk system. Then, we calculate various observational signatures {Tyoi, Lyoi/Lsmrm and 
infrared colors) as a function of time. 

As defined by the bolometric temperature criterion, the Class 0 stage should be very short, while 
the Class I stage persists for much of the protostar’s early life. We present physical distinctions among 
the classes of forming stars and calculate the observational signatures for these classes. Finally, we 
present models of infrared color-magnitude diagrams, as observed by the Spitzer Space Telescope, 
that should be strong discriminators in determining the stage of evolution for a protostar. 

Subject headings: stars: formation, low-mass 


1. INTRODUCTION 

In 1977, Shu presented seminal work in the theory of 
low-mass, isolated star formation HIIllS hereafter 
Shu77). He presented the idea that stars could form 
from inside-out collapse. This model is still important 
today because it is simple, yet it predicts so many ob¬ 
servables in the process of star formation. It prescribes 
the evolution of inflow, the velocity structure of the en¬ 
velope, and the particular shape of the envelope’s den¬ 
sity distribution. iMotte V Andrei ll2flnifl found that the 
inside-out collapse model fit millimeter observations of 
protostars in Taurus and various Bok globules. How¬ 
ever, these authors also found Class 0 sources in Perseus, 
a less quiescent region of star formation, that had cen¬ 
tral densities and accretion rates that were too high to 
be accounted for by the Shu77 model. Molecular line 
observations have also been used to compare the pre¬ 
dicted densities and v elocities with the actual condi tions 
in star-forming co res. iHogerheiide fc Sandelll ||2000D and 
iZhou et all 1)19931) presented evidence that the envelopes 
of some protostars are undergoing inside-out collapse. 
Others have presen ted evidence against inside-out col¬ 
lapse. For example, iTafalla et al.l lll998D suggested that 
L1544, a starless core, exhibits infall motions across the 
entire core. Nonetheless, it is still not clear whether the 
Shu77 inside-out collapse scenario or its variants can be 
ruled out. 

Therefore, we have calculated the observational signa¬ 
tures of a star, which is forming through inside-out col¬ 
lapse, so that astronomers can test this theory in a well- 
defined way. The most relevant aspects of the Shu77 
inside-out collapse model are a constant accretion rate 
of material from the envelope onto the star-bdisk system 
and an envelope density that is initially described by a 
singular, isothermal sphere (SIS). 

The constant accretion rate in this model has given rise 
to the so-called “luminosity problem”: the luminosities 
that result when material accretes onto a central object 
with a small radius and increasing mass exceed those seen 


for m ost young, low-mass stars l|Keuvon fc HartmannI 
dig). The presence of a disk can help by increasing the 
accretion radius and acting as a reservoir where matter 
is stored and then episodically accreted. However, a disk 
does not completely eliminate the “luminosity problem” 
in these models. 

A constant accretion rate in star-forming cores should 
be evident when comparing populations in the different 
stages of star formation. The transition from Class 0 
to Class I object is thought to occur when the mass of 
the star and disk is equal to the envelope mass. There¬ 
fore, we should observe equal numbers of Class 0 and 
I objects if they fo r m wi th constant accretion rates. 
lAndre fc Montmerlel ill 9941 found about a 10:1 ratio for 
Class I and 0 sources in Ophiuch us, suggesting the C lass 
0 stage is very short. However. iVisser et al.l ll2002D . in 
an unbiased survey of dark clouds, found a 1:1 ratio. 
These authors suggested that Ophiuchus has experienced 
a burst of star formation in the past, resulting in the 10:1 
ratio. Future and ongoing surveys of nearby star-forming 
regions will certainly offer more information regarding 
the relative populations of Class 0 and Class I sources 
ijEvans et af .11200^ iBenia.min et af.lf2003r . 

The initial density configuration for the Shu77 model 
is that of a SIS, which has n(r) oc r~^. One¬ 
dimensional modelling of the submillimeter emission 
from starless cores has shown that their density dis¬ 
tributions are well-fitted by Bonnor-Ebert spheres, but 
power-law density profiles are not conclusively rul ed out 
i|Ward-ThomDSon et al.l 119^ lEvans et alJ I2001D . In¬ 
deed, some starless cores have a density structure that 
seems to be approaching the SIS. For example, three- 
dimensional modelling of L1544 show s a power-law den - 
sity distribution with n(r) oc r' ^ llDotv et al.l 1200, ’ll) . 
Other studies of the density distribution in more evolved 
star-form ing cores were unab le to rule out the Shu77 pre¬ 
dictions llYoung et al.ll2003D: in fact, one-third of the 
sources in lYoung et al.l ll200^ were well-fitted by the 
inside-out collapse model. For this paper, we assume that 
starless cores must evolve into an SIS. At this point, col- 
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lapse ensues, and we begin modeling the evolu tion of the 
protos tar. This work has been expanded bv IT ;ee et a 1 .1 
ll2004^ to model the chemical evolution of these star¬ 
forming cores. 

Other authors ha ve predicted t he ob served signatures 
of a forming star. iMvers et al.l l|1998ll . hereafter M98, 
developed a framework on which to calculate various sig¬ 
natures of star formation. Our work has been prompted 
by their efforts. However, the methods and models em¬ 
ployed in this work are quite different from those of M98; 
as a result, our conclusions are different. We discuss 
these differences in this paper. 

Because these models are 1-dimensional, we do 
not consider the role of outflows, flattened envelopes, 
or asymmetric disks; therefore, we probably under- 
estimate the amoun t of short-wavelength radiation. 
iWhitnev et al.1 ll200,‘lL hereafter W03) have included some 
of these complications. However, W03 did not create a 
consistent evolutionary model but, instead, considered 
“typical” protostellar objects of different observational 
classes. Fortunately, these authors were able to explore 
the impact of 3-dimensional effects, which we are unable 
to model. In this way, we consider this effort to be com¬ 
plementary to the work of W03 and compare our results 
to theirs. 

The advent of large surveys such as 2M ASS and the 
Spitz e r Space Telescope’s Le gacy programs llEvans et al.l 
120031 iBeniamin et alJ 1200311 provide vast sets of data 
through which theories of star formation may be vig¬ 
orously tested. In this work, we hope to provide some 
tangible means to study the validity of the inside-out 
collapse model in the role of forming stars. 

2. THE MODEL 

First, we define the framework through which we have 
created this evolutionary sequence. In this section, we 
discuss what has been assumed for the dust opacity, in¬ 
terstellar radiation field, envelope structure and dynam¬ 
ics, and the evolution of the star and disk components of 
the system. 

2.1. Interstellar Radiation Field & Dust Properties 

lEvans et alJ ll2001ll showed that, for starless cores, 
the interstellar radiation field (ISRF) significantly af¬ 
fects both the total observed luminosity and the shape of 
the observed submillimeter intensity profile; even objects 
with a luminous, internal source might attribute mos t 
of their luininosity to the ISRF llYoung et al.l 12004^ . 
If vans et al.1 ll200lll scaled th e ISRF by a c o nstaii t, but 
we have used the opacity of iDraine fc Led lll984fl dust 
to attenuate the ISR F with Ay = 0.5 (see Figure 8 in 
lYoung et alJ l|2004h(B . This method simulates the effects 
of low density material in the environs of a star-forming 
core. _ _ _ _ 

These authors iFvans et in 1201)11 l^irlev et alJ 12001 
lYoiing et alJl200^ also found that the multi-wavelength 
observations were best matched by using the opac ities 
of the dust modeled bv lOssenkopf fc Hennind il994D . In 
particular, they concluded that “OH5” dust, fouiid in th e 
fifth column of Table I in lOssenkopf fc Hennind llI994D , 
was optimal for star-forming cores. Unfortunately, the 
modeled data for OH5 dust do es not include wave¬ 
length s shortward of 1.25 /rm. lOssenkopf fc Hennind 
calculated only the values for the dust opac¬ 


ity (k) and not the scattering and absorption cross- 
sections (cTabs and a scat) as needed b y DUSTY, the ra¬ 
diative transfer progra m we have used lllvezic et al 
llvezic fc Elitzuil l997t). T herefore, we have obtained data 
from iPoiiadTet al .1 (|1994l) . which includes the scattering 
and absorption cross-sections for wavelengths as short as 
0.091 ^m. In Figure Q] we show t he opacities for OH5 
dust and the opacity calculated bv iPollack et al.1 llI994l) 
for dust grains with a radius of O.I ^m at a temperature 
of 10 K (hereafter, PI dust); we have assumed a gas-to- 
dust ratio of 100 and give the opacity of the gas in this 
figure. At short wavelengths, these two types of opaci¬ 
ties are in fairly good agreement; unless r is low, how¬ 
ever, the short-wavelength opacity is not relevant. We 
used the opacity given for OH5 and a scat for the PI dust 
to calculate the absorption coefficient for the OH5 dust. 
Furth er, we used the albedo values given bv iPendletonI 
iTToH Figure 4b) to apportion the opacity due to scatter¬ 
ing and absorption from the 3 /rm ice feature. Finally, 
we have extrapolated the cross-sections out to 3.6 cm, 
as required by DUSTY. For tJscat, we extrapolate by a 
A“'* power-law as expected for Rayleigh scattering. We 
fit the last several data points of the OH5 absorption 
coefficients to determine the ® power-law used to ex¬ 
trapolate CTohs out to A = 3.6 cm. We show the scattering 
and absorption coefficients in Figure ^ 


2.2. Envelope 

For the density structure in the envelo pe, we adopt 
the inside-out collapse scenario llShiill1977l) . This model 
begins with an SIS with a density distribution that is 
proportional to r“^. Through some perturbation, col¬ 
lapse begins inside the cloud and proceeds outward. As 
collapse ensues, the cloud’s density distribution can be 
approximately described by a broken power law: the in¬ 
ner collapsing portion, n oc (indicative of freefall), 

and the outer static envelope, n oc r“^. However, there 
is a transition region just within the infall radius where 
the density profile is significantly flatter than the 
power law. Therefore, we use th e actual solutions to 
Equations II and 12 in iShiJ (IT 97 I . 

When the infall radius exceeds the outer radius, the 
Shu77 solution is no longer valid. Therefore, we adopt 
a density profile with n oc r~^/‘^ and let the mass of the 
envelope (and, hence, the fiducial density) decrease as 
mass is accreted onto the protostar and disk. 

The total amount of mass is constrained by the ef¬ 
fective sound speed (cs) and the envelope’s outer radius 
(co). The models presented herein all begin with cores 
whose initial masses are different. We calculate this mass 
from the following expression: 


M. 


t=o 



( 1 ) 


where G is the gravitational constant, and Cs is the ef¬ 
fective sound speed, 


kT I 2 

A Ty^turb ) 

ptmu 2 J 


( 2 ) 


where k is Boltzmann’s constant, T is the isothermal 
temperature, /r is the mean molecular mass (/r = 2.29), 
and niH is the mass of the hydrogen atom {mn = 
1.6733 X I0“^^g) and Vturb is the turbulent velocity (I/e 































































Formation of Protostars 


3 


Doppler width). We choose T = 10 K and the value for 
vturb such that the turbulent contribution to the sound 
speed is equal to the thermal component (vturb = 0.268 
km s“^); Cs = 0.268 km s“^. 

We calculate the total envelope mass as follows, 

To 

^^env fJjTflfj / 4TTr^n{r)dr, (3) 

Vi 

where and Tq are the inner and outer radii of the enve¬ 
lope. In this paper, radii pertaining to the envelope are 
denoted by the lower-case “r” (n, rmf, and Tq) while the 
radii of the star and disk are denoted by the upper-case 
“i?” (i?*, i?„ Rd). 

Even though the Shu77 model has no specific mass 
scale, we assume that the model is applicable to the for¬ 
mation of stars with different masses. In order to define 
the mass of the core, we truncate the outer radius of the 
envelope, using Equ ation HI Such trun cated envelopes 
are not unheard of; iMotte et al.l lll998ll found evidence 
for truncated outer radii for cores in Ophiuchus. In this 
paper, we consider cores with three initial masses: 0.3, 
I.O, and 3.0 Mq. With our assumed sound speed, these 
masses correspond to outer radii of 1767, 5889, and 17667 
AU. We end our modeling when all envelope material has 
accreted onto the star-disk system. In these three scenar¬ 
ios, this event occurs at 62500, 210000, and 625000 years, 
respectively. The time for collapse varies significantly 
among the three models, but this variation is inherent 
within the inside-out collapse model, assuming similar 
initial conditions. With constant and identical accretion 
rates, lower-mass objects simply form more quickly than 
higher-mass objects. We show the mass evolution for 
each of these cases in Figure |21 

It is not clear whether this model is realistic for cores 
with masses as low as 0.3 M©. The Jean’s mass for a 
core with T = 10 K and density of 10® cm“® is ~0.6 M© 
(Mj = 18M©r^®n“® ®). However, if the cloud is cooler, 
it could be unstable to collapse (e.g., if T = 5 K and 
n = 10®, Mj = 0.2 Mq). In fact, much smaller mass 
cores can be created thro ugh turbulent fragmentation 
llBovd fc Whitwortl]l200,5ll . so the 0.3 M© is probably 
not so absurd. Considering that some evolved protostar s 
are thought to be substellar ilWhite fc Hillenbrandl200Jl . 
we must consider how these objects are formed. 

Defining the envelope’s inner radius is an issue. One 
choice is that the envelope’s inner radius could be equal 
to the outer radius of the disk. However, this disk radius, 
which is defined as the centrifugal radius fSection [2. 3. HI . 
is very small at early times. With a small inner radius, 
the density in the inner region is unrealistically large 
(e.g., n ~ 10^® cm“®). These dense regions cause the 
opacity to become very high. Further, a rotating enve- 
lope becomes flattene d and aspherical at these small radii 
llTerebev et al.llT9^ . so a spherical model is not appro¬ 
priate to these regions. Therefore, we set a maximum 
value for r(100 fim) and calculate the inner radius of the 
SIS (at t = 0) so that ri,(100/rm) is equal to T^ax- 

'^max G27r 11”^ 

Ti = -5— -f — (4) 

(Y* 

' o 

We choose r^a x = 10 and discuss the impact of varying 
this in section 10 The envelope’s inner radius follows 


this formula until it is exceeded by the disk radius. For 
the three cores, the envelope’s i nner r adius is 100 AU at 
the end of collapse (see Section 1^.3.111 . 

We define the infall rate for the case of a non-magnetic, 
centrally pea ked e nvelop e density distribution as de¬ 
scribed in the IShul l|1977li scenario. The rate of constant 
accretion is calculated as follows: 

M = mo^, (5) 

where ruo is a dimensionless constant of order unity. 
Since = 0.268 km s“^, the accretion rate, M, is 
4.8 X 10“® M© yr“^. 

We adopt the prescription of lAdams fc Shul lll986H for 
the infall rate onto the disk and star. These authors 
assume that all mass is accreted onto either the disk or 
star such that M = IVR + Mo where M* is the accretion 
rate of envelope material directly onto the star and Mo 
is the accretion rate onto the disk. These values are 
calculated as follows: 

M^ = M[l-{l-u^f/% (6) 

M^3 = M(l-M*)l/^ (7) 

where rt* is the ratio of the star and disk radii, R^^/Ro- 
In almost all cases, m* becomes very small in a short time 
and, hence, M* also approaches zero. However, material 
also accretes from the disk onto t he star; this process is 
not included in these equations. lAdams & Shul lll986ll 
defined an efficiency factor— rjo-, the fraction of material 
in the disk that will accrete onto the star—so that the 
star could gain mass and be a source of accretion lumi¬ 
nosity. We discuss the implementation of this accretion 
process in the next section. 

2.3. Disk 

Evidence in the form of near- and mid-infrared 
jPadget^_et_ayil999|) and millimeter llMundv et al.lll99^ 
IKitamiira et alJl2002li observations of disks surrounding 
stellar and substellar objects has recently become more 
convincing. Also, the inclusion of a disk in various mod¬ 
els has a significant effect on the interpretations of those 
models. Therefore, to not include a disk in this evolu¬ 
tionary scheme is wholly unrealistic. 

We adopt the disk model developed, in theory, 
JArhuns^hm fc Ladal (1198811 . and, in practice, by 
IButner e^LLTil^^r The density distribution for the 
dust and gas (assuming homogeneous mixing) is defined 
as 

E(A) = So(^) ^ (8) 

where So is the surface density (in gm/cm^) at i?/, a 
fiducial radius. We choose p = 1.5 in accordance with 
the density structure for ve rtical hydrostatic equilibrium 
llChiang fc G!oldreichlll997^ . The mass of the disk, given 
this power law distribution is given by: 

y 27rEoA? „ „ 

Mo = 2-KY.RdR = < 2, 

J ^ P 

Ri 

(9) 

where i?i, Ro, and Rf are the inner, outer, and fiducial 
radii of the disk. 
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2.3.1. Radius of the Disk 

The inner radius of the disk is the dust destruction 
radius defined as: 


Ri = 


I L, 


( 10 ) 


where we define the dust destruction temperature, 
Tdust = 2000 K, and L* is the luminosity of the star. 

For the disk outer radius, Rd, we adopt the centrifugal 
radius that evolves with time as follows ijTerebev et al.l 

dH, 

3 

Roit) = (11) 

Id 

where t is the time and flo is the angular velocity of 
the cloud prior to collapse; other variables are as already 
defined. We set Oo so that, at the end of the Class I stage, 
the disk radius is 100 AU. These angular velocities are 
1 X 10-14, 5.5 X 10-14, and 3.4 x IQ-i^ s-i for the 3, 1 , 
and 0.3 Mq models, respectively. iGoodman et all l|1993ll 
found a range for Go from 9.7 x lO-i® to 1.3 x lO-i^ s-i. 
The upper end of this range is about one-third of Go for 
the 0.3 Mf 7 ^ model, bu t the least massive of the cores in 
IGoodman et al.1 ill 99311 was 0.6 Mq. We assume that less 
massive and smaller cores have higher angular velocities. 


2.3.2. Mass of the Disk 

We evolve the mass of the disk via the expression given 
bylAd ams fc Shill lIT^ : 

Md = (1 — rjjy) f Modt = (1 — r]D)M.T>M, ( 12 ) 

■Jto 

where to is the time when = 1 (i.e., Rp = R ^). We 
assume to = 0. Further, lAdams fc ShiJ 11198611 defines 
A4x> as follows: 

1 

Mv = 7:u* / (1 — u*)i/^M-4/^du. (13) 

We evaluate this exp ression numerically as suggested by 
lAdams fc ShiJ lll986ll . Finally, M is the total mass ac¬ 
creted (i.e. M = Mt). lAdams fc ShiJ l|1986ll give t]d as 
a free parameter; rjo is the fraction of material accreted 
onto the disk that will eventually accrete onto the star. 
To determine a value for rji), we assume that the ratio 
of star to disk mass must be ~ 1/4, which is in accord 
with theoretical work lfil l2nfll . In Figured we show 
this ratio for several values of tjd- We choose ijd = 0.75 
so that Mo/Mt, ~ 1/4 for the 1 Mq model. We apply 
the same criterion for the 0.3 and 3 Mq models and set 
rjjy = 0.7 and rju = 0.9, respectively. In sectional we 
explore the effects of allowing it to vary. 

2.3.3. Luminosity of the Disk 

The temperature distribution for the disk is defined by 
the following, 

T(A) = ro(^) \k), (14) 

where To is the temperature at the fiducial radius, Rf. 
We set q = 0.5, a temperature distribution that decreases 
more slowly with radius than expected for a fiat disk, to 


simulate flaring and disk a ccretion llButner et ahllTo^ 
iKenvon fc Hartma,nnllT98^ . 

The lu minosity of the disk h as several components as 
given by lAdams fc (dig). First, envelope mate¬ 
rial falling onto the disk will act as a source of lumi¬ 
nosity. Then, there is a source of “mixing luminosity” 
that arises, basically, from the mixing of newly accreted 
material and that material already in orbit around the 
star. Finally, lAdams fc ShiJ lll986^ assumed that some 
fraction (ry_D) of the disk material will frictionally dis¬ 
sipate its_j;emaimM orbit al energy and fall onto the 
star. lAdams fc Shul lll986jl give the expression for L^cc 
in equation 33b of their work; we use this equation for 
Race’ which includes all three of these components, in the 
evolutionary model. As the disk radius grows larger with 
time and u^, 0 , their exp ression for L^cc simplifies to 

llAdams. Lada fc ShiJll98^ 

^acc ~ (Ib) 

where To = ■ 


2.4. Star 

For the star, we define several parameters: mass (M*), 
luminosity (T*), radius (i?*), and effective temperature 
(Toff). Each of these quantities evolve over time, and 
they are generally interdependent. 

2.4.1. Mass of the Star 

For the stellar mass, we s ubtract the disk mass from 
the total accreted mass as in lAdams fc Shd lll986ll : 

M* = M - Md = [1 - (1 - riD)Mv]M. (16) 

All variables are as previously defined. Included in this 
equation are two means by which the star gains mass. 
Material accretes directly from the envelope onto the 
star until the disk grows in size and accretes most of 
the infalling envelope mass. Then, the star mostly gains 
material that has accreted onto and is processed through 
the disk. 


2.4.2. Luminosity of the Star 

The luminosity of the star has two components: that 
arising from accretion {L*oo), the dominant source of 
luminosity at early times, and the luminosity due to 
gravitational contraction and deuterium burning (Lphot)- 
S imply, T^, = T* -|- L^ hot ■ 

lAdams fc ShiJ lll986ll calculate the accretion luminos¬ 
ity from material accreting onto the star. Their calcu¬ 
lations include the luminosity from material that falls 
onto the star and the energ y released due to diff eren- 
tial rotation of the protostar. lAdams fc Shul l)1986ll give 
an expression for T*^,, in equation 33a of their work; 
we use this prescription for L*^^. As the disk radius 
grows larger with tiine an d u* —> 0 , this simplifies to 
llAdams. Lada fc 

Kcc^^vWLo. (17) 

lAdams fc 123(11983) define rj^, as an “efficiency factor” 
that dictates how t he star dissipates the energy due to 
differential rotation. lAdams. Lada fc ShiJ lll987ll consid¬ 
ered this value to be a free parameter but chose 77 * = 0.5 
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for their standard model, which we use as well. In sec¬ 
tional we discuss the effects of v pying _ 

Lnh nt was calculated by iD’Antona fc Mazzitellil 
<11994 hereafter DM). We have made linear fits to 
their pre-main sequen ce tracks with opacities from 
IMexa.nder et al.l l|1 989fl (Tables 1 and 5 of DM). First, 
we have fit a power-law in the luminosity-time plane for 
each stellar mass given by DM. For masses less than 0.2 
Mq, where a single power-law is not appropriate, we 
have fit two-piece power-laws to DM’s data. For times 
earlier than those covered by DM’s tracks, we have 


assumed a power-law expression: 


-'phot 


_ j phot 

— i-/o 



where to is the earliest time in the calculations by DM, 
and Lo^°* is the luminosity of the pre-main sequence 
star at time to- This equation is ad hoc and meant to 
smoothly bridge the transition from where there is no 
data for Ijptot to the point where DM’s evolutionary 
tracks begin. Finally, to obtain the appropriate value 
for Lphot, we linearly interpolate between masses for a 
given time in the star’s evolution. 

As noted in M98, the beginning of infall and accretion 
luminosity is not the same time as that for the onset of 
the luminosity represented by Lphot- We adopt, as M98 
did, a difference in these two timescales of 10® years such 
that, for Lphotit), we take t = tphot + W^ yr for particular 
values of tphot as given by DM; th i s assu mption is based 
on the theoretical work of iStahle'rl <1198,8(1 . After collapse 
begins, the forming star must wait 10® years before the 
luminosity due to contraction and deuterium burning (as 
described by DM) will begin. Further, the luminosity 
from DM’s models does not become significant until ~ 
7 X 10"* years. Therefore, Lphot is not relevant until t ~ 
1.7 X 10® years, which is greater than the time required 
for the 0.3 Mq to collapse and only 4 x 10"* years more 
than the collapse time for the 1 Mq core. 

With this final term, the total luminosity of the pro- 
tostellar system is now made of three components: L*^^, 
the luminosity due to accretion and differential rotation 
of the protostar; Lphot, the luminosity arising from grav¬ 
itational contraction and deuterium burning in the pro¬ 
tostar, and luminosity from accretion onto the 

disk and dissipation of the orbital energy within the disk. 

In addition, there is the luminosity that results from 
the ISRF, Ljsrf- We calculated Ljsrf by illuminating 
a core that has no internal source. However, the core 
does have an evolving density distribution identical to 
the three mass scenarios presented herein. Then, we cal¬ 
culate the luminosity that results from the dust grains, 
which are heated externally. At early times, the exter¬ 
nal radiation field contributes most of the luminosity. As 
the envelope mass decreases and the internal source lu¬ 
minosity grows, Ljsrf becomes insignificant. We plot 
the evolution of Ljsrf in Figure 01 

In conclusion, the total luminosity is given as 


Ltot — -^acc + Lphot + -b^cc + Ljsrf- 


(18) 


For Tgff, we use the Stefan-Boltzmann Law, 


- {a4FRl) 


1/4 


(19) 


where L* = L*^^ -|- Lphot- At early times, the effective 
temperature is very low, ~ 100 K because the radius of 


the first hydrostatic core fSection l2.4.3ll is ~5 AU. When 
the stellar radius approaches 2-5 Rq, Tg// becomes more 
stellar-like (~ 3000 K). 


2.4.3. Radius of the Star - Simulating the FHC 

We allow the rad i us of the star to evolve as suggested 
bv iPalla fc Stable"^ <ll991i see Figure 1 in their paper). 
In their calculations, they find that the radius of the star 
rises to about 2-5 solar radii; this result is in accord with 
the assumption of a constant radius at 3 Rq by M98. 
However, the time at which to apply these calculations 
is not so clear. 

In the evolution of a protostar, the early years are 
occupied by the first hydrostatic core. While not yet 
clearly ob served, the first hy d rostatic core has bee n 
predicted <IBoss fc Yorkel 1199,4 IMasunaga et all 1199811 . 
iBoss fc Yorkel <1199,41 concluded, based on some simple 
arguments, that the lifetime of this stage should be 
short, only about 20,000 years. Further. IMasunaga et alJ 
<1199811 have determined that the average radius of this 
core should be about 5 AU. The transition between 
this very large core and the smaller core described by 
the calculations of <IPa,11afc StahleHir^ is not well- 
understood. We have assumed the stellar radius evolves 
as shown in Figure 0 In the beginning, the radius of 
the first hydrostatic core is 5 AU. At t=20,000 years, we 
allow the r adius to decrease from 5 AU to the radius cal¬ 
culated bv IPalla fc Stable^ <(199111 . This transition lasts 
100 years and is described, in our model, by: 


R^AU) = 5 1 


t - 20000' 


0.5 


100 J 

20,000 < t < 20,100 


R 


PS 


( 20 ) 


where t is the time in y ears and Rf^ is the v alue for 
the radius calculated bv IPalla fc Stahle'3 <1199111 . This 
equation is somewhat ad hoc and simply used to simulate 
the transition between the large radius as predicted by 
IMasunaga et al.l <1199811 and th e much smaller r adius of 
the actual star as predicted bv IPalla fc Stahleil <(1991(1 . 

There are consequences for including this large radius 
at early times. Because the centrifugal radius is very 
small and, hence, the disk has not formed, the luminos¬ 
ity at these early times is derived wholly from spherical 
accretion onto the central source. If this central source is 
small, the accretion luminosity can be very high. In Fig¬ 
ure El we show the evolution of the accretion luminosity 
for two scenarios. In one case, we have included the FHC. 
We also plot L*^^ when there is no F HC. In this case, 
the st ellar radius evolves via the data of IPalla fc Stahlen 
<IT^ (i.e., between 2 and 5 Rq). Without a FHC, 
the accretion luminosity rises quickly because there is no 
disk and material is accreting directly onto the star. At 
about 2 X 10^ years, the centrifugal radius has increased 
so that a disk can form. Then, the accreting material 
is processed by the disk causing which arises from 
accretion onto the star, to decrease. 

In summary, we let the mass of the st ar and the ac¬ 
cretion luminosity evolve as defined by lAdams fc ShiJ 
EMl. the rad i us of the star change as predicted by 
IPalla fc Stahled <(1991(1 (except at early times), and the 
luminosity due to deuterium burning and gravitational 
contraction of the PMS star evolve as calculated by DM. 
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The effective temperature, given these other factors, is 
defined by the Stefan-Boltzmann Law. 

3. SIGNATURES 

In this section, we discuss the various observational 
signatures in the evolution of protostellar systems. 

We calculate t he bolometric tempera ture by the pre¬ 
scription given in IM vers fc LaddI i 199,311 . 

Tboi = [C(4)/4C(5)]hi?/fc = 1.25 X 10-“P, (21) 

where C(^) is the Riemann zeta function of argument 
TO, h is Planck’s constant, k is Boltzmann’s constant, 
and the mean frequency, P, is the ratio of the first and 
zeroth frequency moments: 

OO 

u = h/Io,I„, = J i^^S.du. (22) 

0 

In addition, we calculate the bolometric luminosity: 

OO 

Lbol = J iirD^S^diy (23) 

0 

where Si, is the flux density, and D is the distance. 
We set D = 140 pc, suitable for nearby star-forming 
regions. We also calculate Lboi/Lgmm for our models 
(lAndre et alJll99.‘^ . Lsmm is found by integrating Equa¬ 
tion ESI from 350 /xm to oo. 

We also calculate the fluxes that would be seen by 
the photometric bands on the Spitzer Space Telescope 
by convolving the modeled spectral energy distribu¬ 
tion (SED) with the bandpasses for the MIPS instru¬ 
ments. These fluxes do not vary substantially from the 
monochromatic fluxes for the central wavelength of each 
bandpass. To convert the fluxes to magnitudes, we use 
these zero-point fluxes for MIPS bands 1-3 (24, 70, and 
160 fx m), respectively: 7.2, 0.8, and 0.17 Jy llYoung et al.1 
Eoo|). We do not convolve the MIPS resolution element 
with the model but assume that all emission is included 
within this beam. 

4. RADIATIVE TRANSFER 

We us e the radiative tra nsfer code, DUSTY, as devel¬ 
oped bv llvezic et al to calculate the temperature 

distribution in the envelope and the emergent SED of the 
star, disk and envelope. In this section, we discuss the 
effect and treatment of scattering by dust grains in these 
calculations. 

DUSTY assumes that scattering from dust grains is 
isotropic. Longward of 10 /xm, the SED is barely af¬ 
fected because the scattering cross-section {ascat) is sig¬ 
nificantly less than the absorption cross-section (uabs, see 
Figure^. Further, the effect is also minimal at wave¬ 
lengths shortward of 10 /xm when the interstellar radia¬ 
tion held is not included. 

Unfortunately, the assumption of isotropic scattering 
causes some problems when the interstellar radiation 
Held is included. In FigureEl we show the SED for a core 
with a mass of 1 Mq (with rioo^m = !)■ In these models, 
the only heating of the dust grains is externally from the 
ISRF. For the solid line, we include the effects of isotropic 
scattering by the dust grains while the dashed line shows 
the SED without scattering. Both SEDs have a peak 


at submillimeter wavelengths as is expected. However, 
the SED with scattering included also has a peak in the 
near-infrared. Of course, we do not observe strong near- 
infrared radiation from starless cores. At short wave¬ 
lengths, these dust grains preferentially forward scatter 
light, so neglecting the anisotropic nature of the scatter¬ 
ing causes this unrealistic Hux in the near-infrared. 

Our options are either a) neglect the effects of scat¬ 
tering in calculating the emergent SED or b) ignore the 
ISRF. The latter is not really feasible because, at early 
times, the ISRF provides the sole source of heating and, 
hence, ignoring it radically affects the temperature pro- 
Hle for the core. Therefore, we opt for the Hrst alternative 
and ignore the effects of scattering. 

5. AN EVOLVING PROTOSTAR 

With these methods and assumptions, we have calcu¬ 
lated the SED for an evolving protostar. In FigureEl we 
show the SED for particular times in the evolution of the 
core that began with a pre-collapse mass of 1 Mq. We 
have also set Tm ax = 10, xy xx = 0.75, and 77 * = 0.5 (seesec- 
tions l2.2ll2.3.2l and l2.4.2ll . The solid line is the emergent 
SED as observed at 140 pc, the distance to Taurus. The 
dashed line represents the star-|-disk SED; this spectrum 
is for the central source and is the input for DUSTY. 
The bars represent the sensitivit y for the Spitzer S pace 
Telescope c2d Legacy program i|Evans et al.ll200.'^ : we 
have increased the 70 /xm sensitivity by a factor of three 
based on in-Hight performance. The asterisks in the sec¬ 
ond frame are IRAS sensitivities. 

We calculate the observational signatures, described 
in Section El for models with different initial conditions 
and whose evolution proceeds in different ways. We use 
different timesteps for the models: At = 1000, 2000, and 
6000 years for the 0.3, 1, and 3 Mq models, respectively. 
These timesteps are each about 1% of the total infall 
time. 

6. FREE PARAMETERS 

In this section, we explore the effects of various param¬ 
eters in these models. We show how the model changes 
when we use different values for rjo, ??*, and Tmax- We 
Hnd that neither rju or have a large effect on the ob¬ 
servational signatures of the SEDs, but the choice for 
Tmax can signihcantly affect the short-wavelength emis¬ 
sion during all stages of evolution. We adopt these val¬ 
ues: r]D = 0.75, 7 /* = 0.5, and Tmax = 10. 

6 . 1 . r]D 

This factor determines what fraction of the disk ma¬ 
terial will dissipate its energy and accrete onto the star. 
It is relevant in the calculation of the disk mass {Mu, 
equation P^. the disk accretion luminosity (L^ .^, equa¬ 
tion^)), the mass of the star (M,, equation and the 
stellar accretion luminosity (Lacci equation EJ- 

We show the effects of changing rju for the 1 Mq core 
in Figure ini Varying rju alters the evolution of all com¬ 
ponents of the luminosity. Of course, this parameter is 
included directly in the equations for the disk and stellar 
accretion luminosities. Because rju determines the stel¬ 
lar mass, it also affects the luminosity due to deuterium 
burning and contraction {Lphot)- With rju = 0.75, 
each of these are higher than in the other two scenar¬ 
ios {rju = 0.25 & 0.5). 
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Finally, in the left panels of Figure cni we show how 
changing rjD affects observational signatures. We plot 
Tboi and Lboi/Lsmm versus time with r]D = 0.25, 0.5, 
and 0.75. Variations in ? 7 £) have almost no effect on the 
evolution of Thai, and Lboi/Lsmm only varies slightly for 
the different values oi r]r)- 

6 . 2 . 77 * 

This factor is a measure of how much of the total lu¬ 
minosity arises from the central star. In practice, it is 
only relevant for the accretion luminosity, as in 

equation czi In Figure CH we plot the evolution of T^cc 
for different values of 77 *. The value that we choose for 
this “efficiency factor” does considerably change the ac¬ 
cretion luminosity. However, as shown by Figure EH the 
value for 77 * has very little effect on the observational sig¬ 
natures. For the standard model in this paper, we adopt 
77 * = 0.5. 

6.3. TjYiax 

As discussed in section o the inner radius of the en¬ 
velope is set so that T(I 0077 m) does not exceed a certain 
value. Previously, some have proposed that th e condi¬ 
tion o f isothermal ity is violated if r('100um')> I ijlja.rsonl 
[ms). However, iMasunaga fc Inutsukal lll999ll showed 
that T(I00/rm) could be significantly greater than 1 while 
the core still maintains isothermality. 

In FigureEl we plot the evolution of ry (100 /rm) with 
Tmax = 1, 5, 10, and 15. The density profile begins with 
a distribution of n{r) oc r“^, but, as collapse begins, 
the inner region drops to a considerably shallower pro¬ 
file. Because the inner density profile changes, r drops 
considerably. In Figure IT^ we show the SED of a 1 Mq 
core at t = 4 X 10^ and 10^ years with r^ax = 1, 10, and 
15. The model with Tmax = 1 shows more emission at 
shorter wavelengths than those with higher Tmax- The 
SED with Tmax = 10 is almost identical, at all times, to 
that with Tmax = 15. 

In Figure [13 there is a discontinuity at ~ 1 x 10 ® 
years. This is the point where the infall radius exceeds 
the outer radius of t he in f alling : envelope. At this point, 
we changed from the lShul (1197711 solution to an envelope 
whose density is described by n(r) oc This is not 

a perfect transition, however, so the “kink” at ~ 1 x 10 ® 
years is an artifact of the model. 

We explore the effects of changing Tmax on the observa¬ 
tional signatures of our evolving protostar. In Figure lT?! 
we show the evolution of Ttoi and Ltoi /Tsmm for the dif¬ 
ferent values of Tmax- The bolometric temperature is 
most affected by varying Tmax] higher values for Tmax 
slow the transition from a Class I to a Class H protostar 
(as defined by Thai)- Simply, less short-wavelength radi¬ 
ation can escape the cloud when the opacity is higher. 
However, for those models with Tmax > 5, the evolu¬ 
tion of Thai begins to converge. Finally, not surprisingly, 
Lboi/Lsmm is unaffected by changing Tmax- We do not 
significantly alter the mass of the envelope, which dic¬ 
tates Lsmrm nor the accretion processes, which are re¬ 
sponsible for Lf,oi- 

The actual value for Tmax is highly uncertain. We do 
not fully understand the density structure in this tran¬ 
sition area between the envelope and disk nor do we in¬ 
clude the geometrical effects of a flattened envelope in 


this region. We assume that Tmax = 10; higher values 
of Tmax have little effect on observed quantities (i.e., the 
SED and its derived signatures). The peak of a 10 K 
blackbody is 350 /im. Assuming k oc A“^ ®, = 1-5 

when Tioo/xm = 10 . 

7. RESULTS 

With these free parameters set to the given values, we 
have run models with different initial masses: 0.3, 1.0, 
and 3.0 M©. Then, we have calculated the various ob¬ 
servational signatures. 

7.1. Tbol QTld Lbol /L smm 

The bolometric temperature and the ratio of bolomet¬ 
ric to submillimeter luminosity have emerged, in the past 
decade, as the two primary methods of classifying pro¬ 
tostars. Because T^oi is a measure of the flux-weighted 
mean frequency of the protostar’s SED, it is highly af¬ 
fected by the emergence of any short-wavelength radia¬ 
tion. The ratio of the bolometric to submillimeter lumi¬ 
nosity, on the other hand, is virtually unaffected by the 
short-wavelength emission. Therefore, Lboi/Lsmm is less 
susceptible to the effects of geometry, which can cause 
more or less NIR radiation to be observed. This ratio 
is a rough measure of the ratio of protostellar mass (in¬ 
cluding the disk and star) to the envelope mass. 

In FigureEl we plot the two signatures as they change 
for the evolving protostar. For the 0.3 M© core, each of 
the evolutionary indicators increase drastically at about 
2.0 X 10^ years, when the FHC contracts. All of the en¬ 
velope material has accreted onto the star and disk by 
6.3 X 10"^ years, shortly after the central star has con¬ 
tracted from the FHC stage. Therefore, as the central 
star becomes hotter and more luminous, there is little 
material left in the envelope. Then, T^oi increases be¬ 
cause more short-wavelength radiation is observed, and 
Lboi /Lsmm increases because Lsmm —> 0 as the envelope 
goes away. If this model is correct, these low-mass stars 
should proceed through the FHC stage and, almost im¬ 
mediately, be seen as Class H objects. 

The 1.0 and 3.0 Mq objects track one another fairly 
well, up to 10 ® years, in the T^oi plot (Figure El despite 
the fact that they form on different timescales. The 1.0 
Mq core requires 2.1 x 10 ® years for all envelope mass to 
accrete while the 3.0 Mq core requires 6.24 x 10® years. 
However, both cores evolve from Class 0 to Class I at 
about 50,000 years, which is about 1/4 and 1/10 of the 
total infall time for the 1.0 and 3.0 Mq cores, respec¬ 
tively. 

This transition from Class 0 to I is partly due to the 
sudden “turning on” of the central source as it contracts 
to ~ 3 Rq, and accretion luminosity becomes relevant. 
There is, however, another reason for this sudden transi¬ 
tion. As shown in Figure IT^ the Shu77 model exhibits a 
drastic decrease in r regardless of the adopted value for 
Tmax- Initially, the envelope is described by n(r) oc r“^, 
but, as collapse ensues, it changes to and r drops. 
Because t is so low, any substantial source of stellar lumi¬ 
nosity will cause observable short-wavelength radiation 
to emerge from the system. 

These details are actually quite important if one uses 
these evolutionary signatures to derive relative lifetimes 
of the various classes as has been done in the past. If 
one calculates the bolometric temperature for a group of 
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protostars (whose SEDs have been well sampled), there 
should be very few Class 0 cores—conservatively, about 
1/10 to 1/4 of the Class 1 population , but most likely 
a much smaller fraction. iVisser et al ] dmil presented 
their efforts to do such a study; they found approximately 
equal numbers of Class 0 and 1 objects. However, their 
data only included the far-infrared (IRAS) and submil¬ 
limeter fluxes. Analysis of the mid- and near-infrared 
data from the Spitzer Legacy and 2MASS surveys along 
with far-infrared and millimeter observations will almost 
certainly produce different results. With more complete 
sampling of the protostars’ SEDs, very few Class 0 cores, 
by the T{,oi criterion, should remain if this picture of evo¬ 
lution is correct. 

Finally, in the plot of Tboi, the 0.3, 1.0, and 3.0 Mq 
data exhibit a “kink” at about 3 x 10'^, 10®, and 3 x 10® 
years, respectively. This is the point where the enve¬ 
lope’s density distribution is described by a power-law, 
n(r) oc instead of the Shu77 solution. The power- 

law distribution has a slightly higher r that causes less 
short-wavelength radiation to be observed. As a result, 
Tboi decreases slightly at this point of transition, but this 
is an artifact of the model. 

The ratio of bolometric to submillimeter luminosity 
seems to be much more consistent in describing the evo¬ 
lution of these protostars. Adopting Lboi/Lsmm = 200 
as the dividing line for Class 0 and I cores, we find that 
the 1 and 3 Mq protostars become Class I objects after 
1.18 X 10® and 3.6 x 10® years, respectively. These times 
correspond to slightly more than 1/2 of the total infall 
time for each core whereas the T^oi criterion showed that 
the cores became Class I after 1/4 and 1/10 of their total 
infall times. 

Finally, in Figure El we show a plot of T},oi and 
Lboi / Lurn.r rj. for the three models. The poiri t s repr esent 
data fronr lYoung et al.1 ll2003ll . lShirlev et "m ll200(Tll . and 
iFroebrictil ll200,’ill .' In general, these models are consistent 
with the data. However, the 11 starless cores in the lower 
left-hand section of this plot show higher Lboi / Ls^nm than 
the model predicts. The definition of Lboi/Lsmm includes 
data longward of 350 ^m, but, for all of these cores, no 
350 /xm data exist. Therefore, the observed Lsmm is 
lower than that which is modeled. Second, there is lit¬ 
tle near- or mid-infrared data available for many of the 
sources represented. With future observations, the bolo¬ 
metric temperature will almost certainly increase. For 
example, we consider the 1 Mq model and calculate Tboi 
by including different fluxes. If we use IRAS fluxes only, 
Tboi = 83 K for the core at 5x10"^ years while Tboi, cal¬ 
culated with just the Spitzer bands, is 92 K. At t = 10® 
years, Tboi is 88 and 151 K as measured with the IRAS 
and Spitzer fluxes, respectively. 

7.2. Mass Ratio 

We can look at classification from a different, more 
physical, perspective. In Figure HTl we plot the ratio of 
the stellar and disk mass to the envelope mass. Physi¬ 
cally, we might consider a protostar to move from Class 
0 to Class I when this ratio is 1 and there are equal 
amounts of mass in the protostellar system and the en¬ 
velope surrounding this star-fdisk. This event occurs at 
t = 3.1 X 10^, 1.05x10®, and 3.1x10® years for the 0.3, 
1.0, and 3.0 Mq cores. Of course, as defined by Shu77, 
this is also the time when the infall radius is equal to the 


outer radius. 

In Figure El we plot Tboi and Lboi/Lsmm as a func¬ 
tion of (M* -I- MD)/Menv In the Tboi plot, the 1 and 3 
Mq cores change from Class 0 to Class I while (M* -|- 
MD)/Mf.nv < 0.5. With the presently defined bound¬ 
aries, the bolometric temperature does not appropriately 
classify the stages of star formation. Only the 0.3 Mq 
core changes from Class 0 to I when (M* -|- Mjj)/Menv ~ 
1 as is appropriate for our understanding of these stages. 

However, with the Lboi/Lsmm criterion, we find that 
the cross-over from Class 0 to I occurs approximately 
when (M* -|- Mo)/Mf.nv = I, which is a more realistic 
view of these evolutionary stages. Indeed, this observa¬ 
tional signature is also favored because it is not largely 
dependent on what is observed at short wavelengths 
wher e geometric effects play a big role llAndre et alJ 

IT^ . 

Therefore, we set some physical divisions for the evo¬ 
lutionary transitions. First, we let the PPC/Class 0 
transition occur when the FHC first collapses at ~ 
2.0 X lO"* years. For all three cores, this occurs when 
Lboi/Lsmm ~ 35. For the Class O/I transition, we let 
(Af* -|- Mjo) /Menv = I) Lboi/Lsmm ~ 175 when this 
criterion is met. Notice that t his value is s lightly less 
than the requirement given bv lAndre et al.1 lll993ll . i.e. 
Lboi/Lsmm = 200. The Class H stage begins when all 
of the envelope material has been accreted. Our models 
are not reliable at these late times, and the Lboi/Lsmm 
signature is not a very good indicator for this stage of 
evolution. Nonetheless, we find that Lboi/Lsmm is ap¬ 
proximately 2000 at this point, but this value is depen¬ 
dent on the adopted model for the disk. 

7.3. BLT Diagrams: A Comparison with M98 

In FigureEl we show a plot of the bolometric luminos¬ 
ity and temperature, which M98 called a BLT diagram. 
The axes are laid out to mimic the Hertzsprung-Russell 
diagram with Tboi increasing right to left. In Figure El 
we have included M98’s models from Figure 7 in their 
paper. Two of the thin lines in Figure IT^ are their mod¬ 
els for forming, 0.5 Mq protostars whose initial envelope 
masses were 1 and 3 Mq. We also show their 0.3 Mq 
model from Figure 9, but we label it here as 1.8 Mq be¬ 
cause this is the mass of the envelope before collapse be¬ 
gins while 0.3 is the mas s of the star at t=oo. We als o 
plot data fromlYo ung et all ll2003ll . lShirlev et 'HI ll2002|) . 
and iShirlev et al .1 11200411 as c rosses: the dots r e presen t 
data from IChen et al nnn and IChen et~ ] ElS). 
Our tracks are markedly different from any of those pre¬ 
sented in M98, so a summary of the differences between 
two models is relevant here. Primarily, our methods dif¬ 
fer in that M98 attempts to create a reasonable model 
that fits the data while we are simply determining the 
observational signatures of the Shu77 model. 

The most significant difference between this work and 
that of M98 is the assumptions for infall evolution. M98 
described the infall and accretion with an exponential 
decay function such that the accretion rates began at 
about 10“® Mq yr“^ and, as the star approached the 
main sequence, finished with 10“® Mq yr“^. We assume 
constant accretion onto the star-|-disk system throughout 
the duration of the life of the envelope as predicted by 
the Shu77 collapse solution. However, in our model, the 
rate of accretion onto the star’s surface does decrease 
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with time as the disk forms and takes a more prominent 
role in processing material from the envelope to the star. 
Also, the modeled evolution is longer for M98 (10® years) 
than in our model (2 — 6 x 10® years). 

These different assumptions about infall have several 
implications. First, in our model, we form a more mas¬ 
sive star from similar initial conditions in less time. For 
our 3 Mq core, the star reaches 0.5 Mq at t = 132,000 
years. On the other hand, M98’s models require, by de¬ 
sign, about 10® years to create this 0.5 Mq star from a 
3 Mq core. Of course, at the end of 10® years, the star 
created by M98 has completed its pre-main-sequence evo¬ 
lution, while our model still requires time to completely 
accrete the disk material. 

This discrepancy in the star’s final mass presents an¬ 
other difference between this work and that of M98. M98 
included a dispersal timescale for the envelope that in¬ 
cluded an assumption of mass loss due to outflow from 
the central protostar. We do not include outflows in any 
way in these models. The lack of outflows in this work 
has two implications: 1) the mass evolution, as depicted 
in Figure[21 is incomplete and 2) we do not consider the 
effects of an evacuated outflow cavity. Fortunately, the 
mass evolution is not considera bly altered by the exclu¬ 
sion of outflows. iCalved lll998ll found that the ratio of 
mass loss rate to mass accretion rate is ~ 0.1. We are 
unable to model the scattered light coming from the out¬ 
flow cavity. Others have, however, and we discuss their 
work in Section O 

M98 also assume the envelope to have a density profile 
with a free-falling structure, n{r) oc . We use the 
Shu77 solution, which has an inner free-falling envelope 
surrounded by a static envelope with n(r) oc r~^. These 
differences in structure of the envelope cause great dis¬ 
parities in the opacity. For example, a core that has 0.8 
Mq of material with n{r) oc creates ri,(100^TO) = 

0.26. A core with the same amount of material but 
described by the Shu77 collapse solution, with an in¬ 
fall radius that is one-half of the outer radius, creates 
Ti,(100^TO) = 0.17. Such a disparity causes large changes 
in Tboi- For example, if we place a 0.7 Lq star with 
Te//=2000 K inside these two cores, the Shu77 core has 
Tboi = 416 K, and the free-falling core has a bolometric 
temperature that is half as high, Tboi = 207 K. 

M98 do not calculate the full SED. Instead, they con¬ 
sider the optically thin and thick limits and calculate 
two moments of the protostellar spectrum: the bolomet¬ 
ric temperature and luminosity. On the other hand, we 
use DUSTY to calculate the full radiative transfer in the 
protostellar system. 

Finally, M98 use a single power-law to describe the 
dust emissivit y, k. In contrast, we use the du st properties 
calculated bv lOssenkoof fc Hennind While the 

dust opacity is aptly described by a power-law at long 
wavelengths, the shorter wavelength opacities are clearly 
not properly represented in the same way (see Figure [Hi. 

Interestingly, the data encompass both models, but 
M98’s models obviously best cover the median range for 
the data. However, there are a substantial number of 
sources that have a lower bolometric luminosity than ei¬ 
ther model allows. Perhaps, these sources are in the qui¬ 
escent stage of episodic accretion so that Lace is very low, 
or these low-luminosity objects could simply have a mass 


less than 0.3 Mq and not be included in these models. 
Also, our models show higher luminosities at later times 
than most of the data, a result of the aforementioned 
“luminosity problem.” Perhaps, accretion does occur 
episodically and only for short times, and we should ex¬ 
pect only a few objects to be in the phase where mate¬ 
rial is accreting onto the star and, hence, have a high 
luminosity. Another explanation, of course, is that the 
assumption of constant accretion throughout the star’s 
evolution is wrong. 

7.4. Infrared Color-Magnitude Diagrams 

In Figure 1201 we show color-magnitude diagrams as 
would be observed with the Multiband Imaging Pho¬ 
tometer for Spitzer (MIPS) on the Spitzer Space Tele¬ 
scope (SST). The three mass sequences (0.3, I.O, and 3.0 
Mq) are represented by the black lines in these figures. 
We show the magnitude at 24 fim ([24]) plotted against 
the [24] — [70] color; on the right of Figure (201 we plot 
[70] and [70] - [160]. 

Also, we show the models of W03 from Figure 8a of 
their work. The colored lines show W03’s calculations 
over varying angles of inclination. The magenta line rep¬ 
resents almost all inclination angles for the Class 0 stage, 
and the magenta triangle is the Class 0 stage as viewed 
pole-on. 

Many differences exist between these models and our 
work. For example, W03 d escribe the e nvelope as a rotat¬ 
ing, freely falling envelope (jm ric hll976ll : W03 also set the 
envelope’s inner radius to be quite small (~ IORq) while 
we use a much larger inner radius. In addition, W03 have 
a number of “common” model parameters, which do not 
change for the different models. These parameters in¬ 
clude the stellar radius, temperature, and mass as well 
as the overall source luminosity. In our model, we allow 
all of these to evolve with time. Finally, W03 included 
the effects of outflow cavities, 2-dimensional disks, and 
varying inclination angles. We are unable to account for 
these things. 

Young objects in the PPC/Class 0 segment of their 
lifespan are very red in our calculations and occupy the 
lower, right-hand section of this plot. However, these 
models are bluer in this plane than the W03 models for 
a Class 0 source. Several factors probably contribute to 
this disparity. Most of the difference arises from the fact 
that W03’s models have a considerably more luminous 
central source than our models. Thus, the envelope ma¬ 
terial is heated more and emits a greater 70 /rm flux. 
Because W03 uses a very small inner radius, the opti¬ 
cal depth is larger than for our models even though they 
may have similarly massive envelopes. Therefore, rela¬ 
tive to the 24 /xm flux, there is a higher 70 /xm flux for 
W03’s models, and their Class 0 stage appears redder 
(and, of course, more luminous). The effect is much less 
pronounced in the [70] — [160] plane because these longer 
wavelength fluxes are less affected by optical depth ef¬ 
fects. 

As with W03’s models, our evolutionary tracks show 
that the objects become bluer in the [24] — [70] color as 
the source becomes more luminous and [24] increases. 
We have marked the PPC/Class 0 and Class O/I tran¬ 
sitions (as determined by Lboi/Lsmm as in section TT^ . 
The Class I/H transition occurs at the end of our mod¬ 
eling sequence when all envelope mass has accreted onto 
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the star and disk. However, aspherical effects are most 
relevant with these later stages, so these models probably 
underestimate the 24 /rm flux for the Class I/II transi¬ 
tion. Transitions for the 0.3 Mq models are shown as 
cyan squares, the red squares represent the transitions 
for the 1 Mq model, and the green squares are for the 
3.0 Mq model. 

The right panel of Figure 1201 shows a simpler behavior 
for the evolution as seen at 70 and 160 fim. As an FHC, 
the protostar appears very red but becomes bluer as the 
core collapses to ^3 Rq and grows in luminosity. The 
160 /rm flux increases over time as the internal luminosity 
continues to grow. However, there is a point where the 
envelope contains so little material that the 160 /im flux 
drops despite the fact that the far-infrared emission of 
the internal star is growing. By the end of these tracks, 
the objects have reached the end of their Class 1 stage; 
they no longer have envelope material but do have an 
optically thick disk. 

8. CONCLUSIONS 

We have presented the results from modeling the evo¬ 
lution of protostars with 3 different masses: 0.3, 1.0, and 
3.0 Mq. The framework for the evolution of these proto¬ 
stars w as taken mostly from the work of I Adams fc Shul 
(|1986(1 but also used the results from several other au¬ 
thors. These e fforts were s i milar and complementary 
to the work of iMvers et alJ lll998ll and iWhitnev et al.l 
(l2f)f)l but employed different methods and theories and, 
hence, got different results. 

We note that the evolution of these modeled protostars 
is significantly affected by the existence and lifetime of 
the firs t hydrostatic core. The work done heretof ore is 
useful ijRoss fc YorkFirr99fil: iMa.sima.ga, et a,1.ll1 998jl . but 
more detailed theoretical work is needed. As we begin 
to obs erve earlier stages of star formation llYoung et al.1 
I2nn4a|l . the role of the FHC must be ascertained. 

We find that the Class 0 stage, when determined by 
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Tbol, should be short-lived, lasting only about 1/10 to 1/4 
of the protostar’s life. This result is somewhat model de¬ 
pendent, but it should also extend to other models of star 
formation. Therefore, the surveys being conducted with 
the Spitzer Space Telescope and 2MASS should reveal a 
small fraction of Class 0 to Class I objects. However, we 
suggest not using the Tboi criterion for classification. 

We find that the bolometric temperature is a poor dis¬ 
criminator for protostars at early evolutionary stages. 
Instead, we suggest using th e Lboi/Lsmm signature pro¬ 
posed bv l Andre et al.1 ill 9931) . Based on physical grounds, 
we suggest these boundaries for classifying protostars: 
Lbol /Lsravn — 35 for PPC/ClaSS 0, Lbol /Lsrava — 175 for 
Class O/I, and Lboi/Lsmm 2000 for the Class I/H tran¬ 
sition. The latter is largely dependent on the adopted 
disk model. Also, Lboi/Lsmm is not relevant after all en¬ 
velope mass has been accreted (beyond the Class I stage) 
since it is a measur e of the stellar to envelope mass ratio 
ll Andre et al.lll99,^ . 

We have presented several observational tools by which 
the inside-out collapse model can be effectively tested: 
infrared color-magnitude diagrams, plots of the bolomet¬ 
ric luminosity and temperature, and a plot of Tboi and 
Lboi/Lsmm- Large surveys will produce these observable 
quantities for hundreds of young protostars and test the 
theory of inside-out collapse. 
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Fig. 1.— The top plot shows the opacity of the gas with OH5 and PI dust (assuming a gas-to-dust ratio of 100). We have used the 
OH5 dust for our models. However, we used the PI opacities for wavelengths shortward of 1 fim. The bottom plot shows the absorption 
and scattering cross-sections for OH5 dust. 
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Fig. 5.— We show the data used for the radius of our evolving protostar. The radius in the first 20,000 years simulates the first 
hydrostatic core phase. The radius for the remaining time is from Figure 1 of Palla Si Stabler (1991). The thin, vertical lines show where 
the evolution ends for the various masses modeled. 
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Fig. 8. — Each panel in the plot has the SED for a core with initial mass of 1 Mq, a temperature of 10 K, Tmax = 10, r}j:) = 0.75, and 
T}^ = 0.5. The age for the model is labeled in each panel. The dotted line is the input disk spectrum, and the solid line is the observed 
SED as calculated by DUSTY. Sensitivities for the Spitzer Cores to Disks Legacy program are shown as horizontal bars; IRAS sensitivities 
are in the second frame as asterisks. The distance is 140 dc. 
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Fig. 9.— Contributing forms of luminosity are shown here for different values of rjiy. The accretion luminosity (Lacc) is from material 
accreting onto the star. Photospheric luminosity (Lphot) is calculated by DM94 and due to contraction of the central star and deuterium 
burning. The disk luminosity (Lj^)) has several components, which are described in section The sum of these luminosities is shown 

by the “Total.” 
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Fig. 11.— For the 1 Mq core, we plot Lacc for various values of r;*. For our models, we choose r;* = 0.5 in accord with|AJamSj_Jjada^^_S^3 
|T983). 
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Fig. 13.— For the 1 Mq core, we show the modeled SED with different values for r-max- The model with r-max = 1 shows greater 
emission at short wavelengths than those with higher Tmax- We choose Tmax = 10 for the models hereafter. 
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Fig. 14. 


For the 1 Mq core, we show T^yoi and Li^oilLsmm with different values for r-max- 
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Fig. 15.— We plot and Li^oilLsmm as they evolve with time for each of the three models. The horizontal lines represent the 
boundaries for classes of protostars; vertical lines mark the point when all envelope material has accreted onto the star+disk system for 
each of the three modeled cores. 
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kol/lw 

Fig. 16.— We p l ot Thni and Lhni /Lsm .m. fo r the three scenar ios that we modeled. The points are data taken from l^ung et ^ <1200311 . 
IShirlev et al.l 120031 , IShirlev et alJ <1200411 , and IFroebrichl 1200511 . 
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Fig. 17.— For the three models, we plot (M* + The thin, vertical lines mark the point when all of the envelope mass has 

accreted onto the star+disk system. The horizontal line shows where the envelope mass is equal to the accreted protostellar mass. 


TABLE 1 

Model Parameters 



Symbol 

Adopted value 

Envelope parameters 

Envelope inner radius 

n 

Equation 1^ 

Envelope outer radius 

To 

Equation rn 

Envelope dust opacity 

Ki/ 

Ossenkoof & Hennine: ^994") 

Disk parameters 

Inner radius 

Ri 

Equation 1^ 

Outer radius 

Ro 

Equation fffl 

Accretion luminosity 

Ld 

Equation 33b. Adams ^ Shu l'1986') 

Temperature power-law 

Td 

Equation 1141 

Stellar parameters 

Accretion luminosity 

J-^acc 

Equation 33a. Adams A; Shu 1198t)) 

Photospheric luminosity 

Lphot 

D’Antona Xt: Mazzitelli 

Effective temperature 


Equation 1191 

Radius 

R* 

Palla Kr. Stabler ('1991') 
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Fig. 18.— Two evolutionary signatures, and L\,oilLsmm^ are plotted against the ratio (M* M£))IMenv The thin horizontal lines 
denote boundaries for the class transitions while vertical line marks where the star+disk mass equals the envelope mass. 
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Fig. 19.— As seen in M98, this is a bolometric luminosity and temperature diagram. The thick lines are our models for three different 
initial masses. The thin lines are models from M98 with the initial envelope masses labeled. The 1.0 and 3.0 Mm mod els a re from figure 
7 in M98. The 1.8 model is from figure 9 in M98. The crosses are from lYoung et all II20Q3I ). IShirlev et alJ i2Q02f l. and lShirlev et all 
OqoI). The dots are from IChen et in 1199511 and IChen et al.l II1997I) . 
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[ 24 ]-[ 70 ] 



Fig. 20.— The black lines in this plot represent the SST MIPS colors as calculated for our models; [24], [70], and [160] are the magnitudes 
observed with MIPS. The colored lines are from Whitney et al. (2003) and extend over varying inclination angles. For the Class 0 models, 
the line represents mos t inc lination angles except for the pole-on scenario, which is shown here as a magenta, filled triangle. Transitions, 
as discussed in Section o for the 1 Mq models are shown as red squares while the green squares represent the transitions for the 3 Mq 
model and cyan for the 0.3 Mq core. 






